Polyhomogeneous expansions close to null 
and spatial infinity 



Juan Antonio Valiente Kroon * 

Max Planck Institute fur Gravitationsphysik, 

Albert Einstein Institut, 

Am Muhlenberg 1, 

14476 Golm bei Potsdam, 

Germany 

Abstract. A study of the linearised gravitational field (spin 2 zero-rest- mass field) on 
a Minkowski background close to spatial infinity is done. To this purpose, a certain 
representation of spatial infinity in which it is depicted as a cylinder is used. A first 
analysis shows that the solutions generically develop a particular type of logarithmic 
divergence at the sets where spatial infinity touches null infinity. A regularity condition 
on the initial data can be deduced from the analysis of some transport equations on the 
cylinder at spatial infinity. It is given in terms of the linearised version of the Cotton 
tensor and symmetrised higher order derivatives, and it ensures that the solutions of the 
transport equations extend analytically to the sets where spatial infinity touches null 
infinity. It is later shown that this regularity condition together with the requirement 
of some particular degree of tangential smoothness ensures logarithm-free expansions 
of the time development of the linearised gravitational field close to spatial and null 
infinities. 

1 Introduction 

In an analysis of the behaviour of the gravitational field close to null infinity 
and spatial infinity has been given. To this end, a new representation of spa- 
tial infinity was introduced. In this representation spatial infinity is depicted 
as a cylinder, as opposed to the standard representation of it as a point (see 
e.g. pl[). This representation has among other things, the following important 
feature: it allows to formulate an initial value problem in which the data and 
the equations are regular; furthermore, spacelike infinity and null infinity have a 
finite representation with their structure and location known a priori. The afore- 
mentioned analysis shows that a certain type of logarithmic divergences arise at 
the sets where spatial infinity "touches" null infinity (which we will denote by I + 
and I~) unless a certain regularity condition is satisfied by the initial data. The 
sources of these logarithmic divergences can be ultimately traced back to the fact 
that some of the evolution equations degenerate at the sets 1^ (the symbol of the 
system looses rank). This precludes the direct application of standard techniques 
of partial differential equations if one wishes to push the solutions of the field 
equations all the way up to null infinity. This kind of degeneracies is a peculiarity 
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of not only the gravitational field, it is also shared by linear massless fields. In 
order to shed some light on the nature and consequences of these degeneracies, 
here we will look at a simpler situation. Namely, the linearised gravitational field 
(spin 2 zero-rest-mass field) propagating on a Minkowski background. The final 
product of our analysis will be isolation of a series of requirements one needs to 
impose in order to obtain logarithmic-free expansions close to null and spatial 
infinity. One expects, in principle, that it is possible to extend this discussion to 
the gravitational field. 



2 Minkowski spacetime close to null and spatial infinity. 

Our point of departure is the usual representation of Minkowski spacetime M in 
Cartesian coordinates: 

9 = nnvdy^dy", 

where rj^, — diag(l, — 1, — 1, —1). We are interested in analysing the geometry 
of this spacetime close to both null and spatial infinities. Intuitively one expects 
this region of spacetime to contain the domain D = {y^y^ < 0}. Let us start by 
considering the inversion in D given by, 

y A y\ x*x x 

This inversion clearly maps D onto itself. Observe how a point that is far from 
the origin of the y^ coordinates seems to be close to the x^ origin. In this sense, 
the x^ coordinates convey the notion of "lying at infinity" . A simple calculation 
shows that: 

rj^dx^dx" . 



" (x*x x ) 

Let now p be the standard radial coordinate associated with the spatial coordi- 
nates x a , a = 1,2, 3. Then x x x\ — (x ) 2 — p 2 . Finally, the introduction of a new 
coordinate r (|r| < 1) via x° = pr yields the Minkowski metric in the form, 

9 = p 4 (1 ^ T 2)2 ( d (P T ) 2 - d P 2 - P 2d(j2 ) ' 

where da 2 is the standard line element of the unit 2-sphcrc in spherical coordi- 
nates. The latter line element suggests the introduction of two different conformal 
factors. Namely, 

~=P 2 (1-T 2 ), 

= p(l-r 2 ). 

By means of these two conformal factors one obtains two different extensions M 3 
and of the original Minkowski spacetime near spatial infinity. The choice of 
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the first conformal factor a yields the standard representation of spatial infinity 
as a point. Consider the conformally rescaled metric: 

9e = Z 2 g = d{p T f - dp 2 - p 2 da 2 , (1) 
= (da; ) 2 - (dx 1 ) 2 - (dx 2 ) 2 - (dx 3 ) 2 . 

Here and all throughout this chapter tilded quantities will always refer to quan- 
tities in the physical (i.e. not conformally rescaled spacetime). It can be shown 
that all radial spacelike geodesies in M map to spacelike curves (not necessarily 
geodesies) in Ms with {p — 0} as endpoint. From (|l|) we can see that the set of 
points for which p = is in fact a point, which we denote by i°. On the other 
hand, consider the rescaled metric: 

ge = ® 2 9 = ^ (d(pr) 2 ~ dp 2 - p 2 da 2 ) . 

It can also be checked that all radial spacclikc geodesies in M map to spacelike 
curves with endpoints at {p = 0}. Now, the metric seems to be singular at p = 0, 
but introducing a new coordinate r = —log(p) one gets: 

g e = dr 2 - 2rdrdr - (1 - r) 2 dr 2 - da 2 . 

Observe that ge\ p= n — dr 2 — da 2 , thus the set {p = 0} has now the topology 
of a cylinder with S 2 as spatial sections. We define / = {|r| < l,p = 0}, 
/+ = {r = l,p= 0}, I- ={t = -1, p = 0}, 1° = {t = 0, p = 0}. 
Let u, v be the null coordinates given by 

u = p(l + t), v — p{l — r). 

The curves given by u — u', v! a constant and fixed angular coordinates are null 
geodesies for which it can be checked that map to outgoing (future oriented) null 
geodesies in ML Similarly, the null geodesies v = v' , with v' constant, and fixed 
angular coordinates map to incoming (past oriented) null geodesies in ML Thus 
the set J? + = {t = 1, p ^ 0} from which the geodesies with u — constant start 
corresponds to future null infinity. Similarly, the set = {r — — 1, p ^ 0} from 
which the geodesies with v = constant emanate corresponds to past null infinity. 
The sets 1^ where spatial infinity "touches" null infinity are in a particular sense 
(to be discussed later) special. Therefore they are considered neither belonging 
to the spatial infinity cylinder I nor to J'^ . Notice that the conformal factor 
vanishes on J^, / and 1^. 

Most of the discussion here presented will be concerned with the extension 
Mg> of Minkowski spacetime produced by the conformal factor = p(l — t 2 ). 
Using the null geodesies discussed in the previous paragraph one can construct an 
adapted null (NP) tetrad such that the real vector I is tangent along the incoming 
null geodesies, and the vector n points along the outgoing null geodesies 0. The 

1 The use of the vectors I and n here is reverse with respect to what it is the standard 
notation The vector I is generally chosen to lie along outgoing null geodesies. 

This is done to agree with the spinorial notation used in (n. 
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Fig. 1. Spacetime close to spatial and null infinities: to the left the standard repre- 
sentation of spatial infinity as a point i°; to the right the representation where spatial 
infinity is enviewed as a cylinder. 



complex vectors m and m are then usually chosen so that they span the tangent 
space of the two dimensional spheres S p ^ T i = {p = p',r = r'}, with p 1 , r' given. 
This construction has the problem that m should vanish at least in a point on 
Sp' tT '. For technical reasons it will be therefore convenient to coordinatise the 
2-spheres S p _ T by means of unitary representations of SU (2, C). This can be done 
as the unit sphere S 2 can be identified with SU(2,C)/U(1) which is obtained 
from S 3 ~ SU(2, C) via the Hopf map. Any real analytic function / on SU (2, C) 
admits an expansion of the form, 

oo m m 
/ = / , / ] f"i,k,jT m fc 7 - , 

m=0 j=0 k=0 

with complex coefficients f m ,k,j- The functions T m *. m = 0,1,2,..., j, k = 
0, . . . , m can be related to the spherical harmonics Yi t „ ||. The function / will 
be said to have spin-weight s if its expansion in terms of the T m functions is 
of the form, 

oo 2q 

f = f i^la q-s- 
q=\s\ fc=0 

Important for later discussion will be the fact that any C°° function <f> = <fi(p, t) 
on Mq" x SU(2, C) such that for fixed p it is of spin-weight s, with s independent 
of p has a normal expansion of the form, 

<t>= J2 J2 Y.^ kT 2 q k q -sP p ^ 

p>\s\ q=\s\ k=0 

with cpp-q } k complex. This last result will be extensively used. For further details 
on the T m k - functions, the reader is referred to ||. 
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3 Linearised gravity in the F-gauge 

We will describe the linearised gravitational by means of a spin-2 zero-rest-mass 
field 4> ab cd satisfying, 

V aa '4>abcd = 0. 

Using the null tetrad described in the previous section, one can recover the {D, 8) 
and (A, 8) equations of the NP formalism (see e.g. (|ll|Jl2| ) 0: 

(1 + r)d T (f> n - pd p 4> n + X + (f> n+ i = (n - 2)<f> n , (2) 
(1 - r)d T (f) n+ i + pd p cj) n+1 + X-(f) n = (n - l)0„+i (3) 

where n — 0, 1, 2, 3. The coefficients <j> n can be shown to have spin-weight 2 — n. 
The operators X + and X_ are complex linear combinations of left invariant 
vector fields on SU (2, C), and can be related to the S and <3 operators of Newman 
& Penrose ||. The NP equations (||) and (||) imply a set of five propagation 
equations, 



(1 + t)cV0 o - pdpfo+X+q)^ -20o, (4) 

d T <f>i + iX_0 o + |X+0 2 = -0 1; (5) 

d T 2 + ±X_0 X + ±X + 3 = 0, (6) 

d r 03 + |*-02 + ±X+0 4 = 3 , (7) 

(l-r)9 T 04 + /o9p04 + ^-03 = 20 4) (8) 
and a set of three constraint equations, 

Td T <j>x - pdpfa + \X+4> 2 - ±X_0 O = 0, (9) 

Td T <t>2 ~ pdpfo + \x + fo - \X^ X = 0, (10) 

Td T <f> 3 - pdpfa + ±X+0 4 - ±X_0 2 = 0. (11) 



This latter set of equations gives rise to equations on the initial hypersurface 
S = {t — 0} which correspond to the constraints coming from the linearised 
Bianchi identities. 

Using the equations (0)-(P|) one obtains a symmetric hyperbolic system |^,^| 
by simply multiplying (H)-(0) by a factor of 2. The resulting system is of the 
form: 

A°d T <P + A 1 d p $ + A + X+$ + A~X-<I> + B$ = 0, (12) 

where A 11 p = 0, 1, 2, 3 and B are 5x5 matrices. It can be more concisely written 
in the form of (A^apd^p + B a = 0. Crucial for our discussion will be to realize 
that the matrix A = diag(l + r, 2,2,2,1 — r) looses rank at r = ±1. This 

2 Here again, in order to agree with reference our notation is reversed with respect 
to the standard one of the NP formalism. Our <f>o, cj>i, 02, 03 and 04 correspond 
respectively to the 04, 03, 02, 0i and 0o of the standard NP notation. 
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degeneracy will be the source of most our problems as this fact precludes the 
direct use of the standard theory of symmetric hyperbolic systems. 
Let Q be the 5x5 matrix with components given by, 

Then the characteristics of the system (|l2|) are the hypersurfaces (j) — k, with k 
a real constant, and such that the scalar field (j> satisfies, 

det Q = 0. 

It can be readily checked that the hypersurfaces u = k and v = k', with u 
and v defined as in the preceeding section are characteristics of the symmetric 
hyperbolic system ([l2]). In particular, the hypersurfaces defined by v = and 
u = with p ^ correspond to and J?~ respectively. Hence, and 
are both characteristic hypersurfaces of the system The evaluation of 
the NP equations (^) on ^ + gives rise to transport equations which enable us 
to calculate the value of the components . . . ,04 on from a knowledge 
of the radiation field 4>q on . Similarly, using equations (|^) one obtains a 
corresponding set of transport equations on/" by means of which it is possible 
to calculate the values of </>o, . . . ,03 from a knowledge of 4 on . 

A further look to the system (|l2|) reveals that at p — the whole system 
reduces to transport equations which enable us to calculate the value of </>„ 
(n = 0, . . . ,4) from their value at the initial hypersurface S. Hence, we say that 
I is a total characteristic of (|l2|). As a consequence, no boundary data can be 
prescribed on /. 

3.1 Initial data for linearised gravity 

That the spin 2 zero-rest-mass field <p a bcd can be used to describe the linearised 
gravitational field can be seen as follows. Consider a family g^(X) of spacetimes 
depending in a C 1 fashion on the parameter A, and such that at A = one 
obtains Minkowski spacetime. Therefore, 

g^u = rjtxu + A/ijuv + C(A 2 ). 

The symmetric rank 2 tensor describes the first order deviation from flatness 
of the metric 5^^. The computation of the curvature to first order in A yields 

s 

Kfj, v \ p = 2V[ /i V|[ l/ /i A ]| p ], 

where V M is the covariant derivative of the flat spacetime. Then the linearised 
field equations take the form: 



K, vX P = 0. 



(13) 
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The tensor Kp V \ p is trace-free and possesses the same symmetries of the Ricmann 
tensor. It can be described in spinorial language by a totally symmetric spinor 
4>abcd satisfying the spin-2 zero-rest-mass field equation: 

V aa '4> abcd = 0. (14) 

Equation ( [li] ) is a sufficient condition for the tensor K^xp to be derivable (lo- 
cally) from some symmetric tensor hp V [ p"3| . If the spinorial field (j> a bcd is set to 
transform as <p a bcd — ^~ 1 4>abcd then the equation (|l4|) is conformally covariant. 

Let us assume that the family of metrics (^(A) arise from the Einstein 
evolution of a corresponding 1-parameter family of asymptotically flat, time 
symmetric initial data given by the 3-dimensional metric, 

It will be assumed that ^g a /3 satisfies the time symmetric vacuum constraint 
equations. Then, the symmetric tensor ^h a p is a solution of the corresponding 
linearised vacuum constraint equations. In the present discussion we will consider 
1-parameter families of 3-metrics ^ 3 'g a f) such that their suitably conformally 
rescaled counterparts ( 3 'g a /3 extend analytically near i, the infinity of the initial 
Cauchy hypersurface. For conceptual reasons it is convenient to distinguish i 
from i°, the spatial infinity of the whole spacetime. It can be proved fij that 
under this assumption the components of the conformally rescaled Weyl spinor 
of the family of metrics <7 M! /(A) in the unphysical spacetime evaluated on the 
initial hypersurface S are of the form, 

oo 

lMA)|;sr= E ^W\s P n , 

p=|2— n| 

where 

P 2q 

q=\2-n\ k=0 

with complex coefficients c n , P -q,k n = 0,1,2,3. Now, upon linearisation one ob- 
tains a similar behaviour for the components of the spinorial field 4> a bcd- Namely, 

oo 

cb n \r=0 = £ r n \r= p n , (15) 
p=\2—n\ 

where 

P 2q 

4>n\r=Q = dn,p;q,kT 2a ' g-2+iV 

q=\2-n\ k=0 
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The components of the Weyl spinor ip n satisfy the constraints coming from 
the Bianchi identities. Hence, the components <p n satisfy automatically their 
linearised version — equations (||)-([ll]) on S = {t = 0}. 

Later considerations will make use of a particular spinorial object: the Cotton 
spinor (sometimes also called Bach spinor). The Cotton spinor B ab cd is the 3- 
dimcnsional analogue of the Weyl spinor || . It locally characterises conformally 
flat 3-metrics in the sense that it vanishes identically if and only if the 3-metric 
is locally conformally flat. One can construct its corresponding linearised version 
bated- Linearisation around Minkowski yields the following expression, 

b abcd = 2D e{a n(t> bcd) c + QD e{a (t) bcd) \ 

where fl = p 2 is the conformal factor of the 3-metric of the Minkowski initial 
data, and D ab its corresponding spinorial covariant derivative. 

Inspired in a similar result by Friedrich we present now the following rather 
technical result. It will be of much use in later discussions. 

Lemma 1. The coefficients dj iP;P: k of the spin-2 zero-rest-mass field satisfy 
the antisymmetry condition 

do iP - p .k\r=o = -dA,p;p,k\ T =Q, P = 0, 1, fc = 0, ...2p. (16) 
Furthermore, 

dj,p; P ,k\T=o = 0, p = 0,...,s*, k = 0,...,2p, j = 0,...,4. 
if and only if 

D (a q b q ■ ■ ■ D aibl b abcd) {i) = 0, q = 0, . . . , s*. 
The proof follows from direct linearisation of theorem 4.1 in 0. 

4 A regularity condition at spatial infinity 

We now proceed to carry out an analysis of the transport equations one obtains 
upon evaluation of the held equations (0)-(^) and (P)-(|ll]). Consider the equa- 
tions and (pl)-(|ri"|). Differentiating them formally with respect to p and 
evaluating at the cylinder at p = one gets: 



{l+T)d T <t? Q +X + <tP x -{p- 2)0o = O, (17) 

9 T ^ + i(X_^ + X + <g) + 1 =O, (18) 

d T f 2 + l{X_c^ + X + ^) =0, (19) 

(l-T)d T c%+X- ( % + (p-2)(% = 0, (21) 
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and 

rd T cft + \{X+<% - X_<f%) - pfi = 0, (22) 

Td T f 2 + \{X+<% - X-<t>\) - V <t>\ = 0, (23) 

rdAl + h( x +<% - x -<&) - P$ = o, (24) 

where we have set 0£ = dp 4> n \p=o- One can expand these coefficients using the 
functions T k l m in the form: 

p 2g 

'Pn= J2 a ^P^ T 2q k q -2+n- ( 25 ) 

q=\2-n\ fc=0 

The coefficients a n ,p;q,k are, in principle, complex functions of r. Using these 
expansions and the equations (18)-(^0|) and (p^)-(p^), one can calculate the 
coefficients ax,p-,q,k> a 2,p-,q.k, a 3,p-,q.k from a knowledge of the ao )P; q,fc and a4 ]P;9j fe 
via a first order algebraic linear system. Some more algebra leads to, 

(1 - T 2 )d 0tP . q .k + (4 + 2(p - l)r)do,p ;g ,fc + (q + p)(q - P + l)ao, P;9 ,fc = 0, (26) 
(1 - T 2 )a 4 ,p; 9 ,fc + (-4 + 2(p - l)T)ct4,p ;9 ,fc + (q + p)(q - p + l)(H,p;q,k = 0, (27) 

for p > 2, 2 < q < p, the overdot denoting differentiation with respect to 
r. The equations (|26| ) and (^) are examples of Jacobi equations. A canonical 
parametrisation for this class of ordinary differential equations is: 

£>(„ a p) a = (1 - T 2 )a" + {/3-a-(a + (3 + 2r)}a' + n(n + a + (5 + I) a = 0. 

(28) 

In our case the parameters are given by: ao — At = —p — 2, 0o = a± = —p + 2, 
n% = p + q, and n 2 = p — q — 1. Regular solutions for equations (|2(|) and fp7|) 
exist for q ^ p, and are given quite concisely in terms of Jacobi polynomials ||14| . 
Namely, 

ao, P; „ fc = C^-_ P g Z 2 r P+2) (r) + C 2 ( V)^ 2 P i-'- p+2) (r), (29) 
a w = D lP ^- p - 2 \r) + £ 2 ( i T 1 )" 2 ^fe 2 '- ? - 2) (T), (30) 



where Ci, C2, -Di and £>2 are constants which can be determined from the initial 
data at r = 0. In the case q = p, the use of some identities of the Jacobi equation 
leads to: 



«w = ( ^7- r 2 (t) p+ ~ ( f ° + Fi f (i +B) ^i- a) p-i ) • (32) 
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The use of Taylor expansions in the integrals shows that they give rise to loga- 
rithmic terms. More precisely, recalling that: 

A* ln(l -r) + p ± +■■■ + —— + A Q 



(l±s)H(l Ts )f+3 ^ / ( 1 _ T ) P ±2 (1-t 

B PT 2 Bi 



B* ln(l + t) 



(l + r)w2 (1 + r)' 



where the ^4's and B's are some constants. From the latter expansions one can 
conclude that the only non-regular coefficients in ao, p; g,fc and a^q^ in the q = p 
case are of logarithmic nature. Analytic solutions arise if and only if ao,p;p,m(0) = 
a4, P ;p,m(0). That is only the case if the constants E\ and F\ in equations ( |3l| ) 
and (|3^) are both zero. Now, using the antisymmetry condition ( [l6| ) of lemma 
[l] one learns that in fact ao,p ; p, m (0) = a4, P;p ,m(0) = 0. It is not hard to see that 
if this is the case, then a^ p;Pim (0) = for j = 0, . . . , 4. Finally, using the second 
part of lemma [l] one can relate this behaviour of the initial data to the vanishing 
of the linearised Cotton tensor and its symmetrised derivatives on i. The results 
of the previous discussion are summarised in the following theorem. 

Theorem 1. The solutions of the transport equations on I corresponding to the 
system and extend analytically to if and only if the regularity 

condition 

D(a s b s ---D ai b 1 b a bcd)(i) = 0, 3 = 0,1,... , (33) 

holds. 

A peculiarity of the logarithms appearing in the solutions of the transport equa- 
tions is that they only occur at the highest spherical harmonics sector at each 
order. This will be of importance later when discussion the asymptotic expan- 
sions of the field <fi a bcd close to J? + . 

The aforediscussed solutions of the transport equations allow us to obtain 
normal expansions of the form, 

<t>n= £ r n P P , (34) 
p>|2-nj 

with the coefficients <\P n given by (^5|). Now, it is of interest to see how the 
form of these normal expansions and the regularity condition (^3|) reflect on the 
structure of the asymptotic expansions close to null infinity. In order to do this, 
one has essentially to reshuffle the normal expansions. So far, our discussion has 
applied to both future and past null infinities. The forthcoming analysis will 
without loss of generality be focused on Nevertheless, a totally analogous 
treatment can be performed for J^~. 

Putting together the results from the analysis of the transport equations 
(p"7|)-([?l|) and (22)-(pi|) one obtains normal expansions for the coefficients <fio 
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and </>4 which are of the form, 

0o=£ (S[0](1 - r) + oS ln(l - r)P[p + 2, 2p] (1 - r)) ff, 

p=2 

04 = ^ T ) + 6 * ln (! - T ) P ^ - 2 ' ^K 1 - I")) P P - 

p=2 

By S[n](a;) it will be understood a generic infinite series in x starting with 
x n . Similarly, P[ni, na](a;) denotes generic polynomials in x of order na and 
whose lowest order term is x ni . The coefficients a*, 6* and those in the se- 
ries/polynomials are given in terms of the functions T m k -. Matters of conver- 
gence of the normal expansions will be addressed in the next section. Now, a 
careful reshuffling in order to obtain expansions in (1 — r) reveal that, 

^ = Y / S[2](p)(l-rr+^2ls[2}(p)+Y^ oyin(l-r) (1 - r)", 

p=0 p>4 \ s=[t/2] J 

P>0 \ s=[t/2] / 

where [t/2] = t/2 if t is even, and [t/2] = (t + l)/2 if i is odd. Notice that the 
component 04 happens to be the most singular one of the field. If the regularity 
condition (p3) is satisfied up to s = s*, then it is not hard to see that, 



J2 Cfc (l - r) fe + d..+i(l - r) s * +1 ln(l - r) + • • ■ , (35) 



^4 = 

fe=0 

where the coefficients Ck are C°° functions of p and the angular variables, and 

2s, +6 

d s ,+l = P S * +3 X! Ds *+ 1 - mT 2s,+6 m s,+5' ( 36 ) 
m=0 

with D m complex constants. Thus, the logarithmic term in the expansion 
has only dependence in the highest spherical harmonic sector possible at this 
order, i.e. q = s* + 3. Finally, if the regularity condition holds also for s — s* + 1 
then d s „+i = 0. 

5 Polyhomogeneous expansions 
5.1 A substraction argument 

So far, nothing has been said about the convergence of the normal expansions 
( |34| ) calculated in the previous section. We will now discuss how this can be 
done. As before, we write # = (<po, <f>\, 4>2, 03, 4>a)- Let, 

N 1 

p=0 1 
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be the N order partial sum. The linear field equations and (p|)-(pd|) are 

such that <£jv itself happens to be a solution of the field equations with truncated 
(to order iV) initial data. The regularity of <1>n will depend on whether the initial 
data satisfies the regularity condition (|3^) to a given order or not. For example, 
if the condition is satisfied up to s = N, then no ln(l — r) term will be present 
in the coefficients <P P , and thus <Pn will be C°°. How could one estimate the 
rest? Recall that the symmetric hyperbolic system derived from equations (Q)- 
(0) breaks down precisely at r = 1. This is rather unfortunate, as we are mainly 
interested in observing the behaviour of the solutions on null infinity. 

Bounds for the rest can be obtained by considering the conformal extension of 
Minkowski, Ms, in which spatial infinity is represented by a point i . The spin 2 
zero-rest-mass field equations given in spherical coordinates are formally singular 
at p — in this representation. Therefore, we resort to Cartesian coordinates 
In this way, the propagation equations read, 



<9o</>o - 


d 3 </)o 4 


-9i0i4 


id 2 <t>i = 


o, 




2c\,0i - 


f di (j>2 


-1- di4> ■ 


f id 2 4>2 - 


- «9 2 o 


= 0, 


2<9 O 02 - 


f <9l</>3 


4- dicfii - 


f i<9 2 </>3 - 


- id 2 (f>i 


= o, 


2<9 O 3 - 


f di 4>4 


4 di4>2 - 


f i<9 2 </>4 - 


- id2<\>2 


= 0, 


do4>i 4 


d 3 ^4 4 


■ <9i0 3 - 


id 2 (f>3 = 


o, 





-1— n . 



while the constraint equations are given by, 

29 3 01 + <9i</> - difo - id 2 (j)a - id2<t>2 = 0, 

29 3 2 4 d\(j)i - d\4> 3 - i<9 2 </>i - id 2 (t)3 = 0, 
2<9 3 03 + difo ~ <9i<^4 ~ id 2 ^>2 ~ id2$± = 0. 

The components of the linearised Weyl tensor rescale as follows: 4> n = p 
n = 0, . . . ,4, where p 2 = (x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 . Thus, the smooth data in 
discussed in §3.1, and such that O = 0(p 2 ), (f>i = 0{p), & = 0(1), </> 3 = 0{p) 
and (f>4 — 0(p 2 ), becomes singular in the Ms picture at i (i.e. p = 0). However, 
if one provides initial data with sufficiently fast decay at i then one can obtain 
energy estimates of the form, 

\\nl* { u)<C\\$\\l, {s ^ (37) 

where || i I ^2 (o) an< ^ II Hl 2 (2;i) denote the L 2 norms over the region 15 and 
the hypersurface E\ respectively. Higher order estimates for the derivatives can 
be similarly obtained. The domain 13 is as shown in figure (Q). Notice how in 
this way one can obtain estimates that go up to ^ + — and through it. In this 
spirit consider, 



X = <P-$ 



N ■ 



3 The symbol 13 is to be read "mho" . 
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T=l 





I 



P 



Fig. 2. The region of integration used to obtain the energy estimate (f37[). The first 
picture shows the domain of integration in the representation where spatial infinity is 
enviewed as a cylinder, and the second in the in the standard representation of spatial 
infinity as a point. Note that £2 which seems to be an extended 2 dimensional surface 
is in fact a point: i°. The surface S3 is spacelike and almost flat, and may reach and 
go through null infinity. 

A close look to the rescaling rules will convince us that: 



It can be readily checked that such a decay implies X\ t= q £ H N ~ A {Ei) 1 where 
H (Si) denotes the fc-th Sobolev space over the domain E\ . Standard theory 
of linear of symmetric hyperbolic systems (see for example [ |l5| ) then guarantees 
the existence of X such that X £ C N ~ 7 (15). Hence, we have estimated the rest 
of the partial sum <&n- If one wants the solution to be of class C k (l3), k 6 N, 
then the previous discussion shows that one needs the regularity condition ( |33| ) 
to hold up to order s = k + 7. 

The disadvantage of this approach is that it is hard to extend for the case of 
quasilinear equations. Therefore, it will be of little use when we eventually try 
to address the full non-linear gravitational field. We need to devise something 
else. 

5.2 An investigation of expansions close to null infinity 

The analysis of the solutions of the transport equations on the the spatial in- 
finity cylinder / show that logarithmic divergences will generically arise at I ± 
for smooth initial data of the sort considered in §3.1 unless the regularity con- 
dition ([33]) is satisfied. Accordingly, if the regularity condition is satisfied, then 
the solutions to the transport equations extend smoothly through / + . This is a 
hint that the evolution of generic smooth initial data will be polyhomogeneous, 
i.e. its expansions will contain ln(l — r) terms. An analysis and a discussion of 
the properties of such spacetimes close to null infinity, including the existence 
of an intriguing set of conserved quantities analogous to the Newman-Penrose 
constants has been given in ]l|, ^6| , p^|JT^ , |20t . The regularity condition given in 
theorem [j] ensures that the solutions of the transport equations (p7|)-(^l[) and 
(p2|)-(p4l) extend analytically to 1^. Now, with the discussion on I + more or 



X| T=0 = 0( P N - 5 ). 
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Fig. 3. The asymptotic characteristic initial value problem close to spatial infinity. 



less settled, one would like to study the effects the degeneracy of the field equa- 
tions on I + has one null infinity. In particular, one would like to know what do 
the logarithmic terms on / + imply on J^ + . How do they propagate? In order 
to investigate this point, we resort to some ideas coming from the asymptotic 
characteristic initial value problem p|,fI|,|i"o[ . Consider a future oriented null hy- 
persurface A/o intersecting at Zq . One can provide characteristic initial data 
in the following way: 04 is prescribed on A/"o, 4>o on the portion of J? + comprised 
between Zq and / + , and (f>2 and <f>i on Zq. Theorems providing existence in 
a neighbourhood 15 C J~(Zq) of Zq for the conformal Einstein equations have 
been given in |J3|j4| Jl 0(1 . An existence theorem for polyhomogeneous Maxwell fields 
has been given in |1°| . Problems with the characteristic approach have mainly 
to do with the ample freedom one has in choosing the initial data on A/"o- As an 
added factor, the theorems by Friedrich and Kannar given only local existence. 
This is somehow a problem, as one would like to take the field 4> a b c d all along 
JP + and evaluate at the set I + , and then compare with the results obtained 
from the transport equations on /. 

Inspired by the aforementioned characteristic approach, we will make use of 
the transport equations obtained on J? + arising from the NP equations ||). In 
this way once the radiation field (f>o is provided, the remaining field components 
can be obtained essentially by a mere integration along the generators of null 
infinity. Expansions containing ln(l— r) terms appeared in the normal expansions 
(f{4|). However, it is of interest to consider a more general type of situation 
in which terms of the form ln fc (l — r) can arise. So, we will assume that the 
components of the spin 2 zero-rest-mass field 4>abcd have asymptotic expansions 
near to J? + of the form, 

^n=EE^ ,i) ( 1 -^^'(l-T), (38) 
i>o j=a 

= 0, . . . , 4 , where the coefficients (f>n contain only p and angular dependence. 
They are assumed to be smooth at least in {r = 1,0 < p > po} for p a 
given non- negative real number. Some calculus rules will be needed in order to 
manipulate the polyhomogeneous expansions (|38|). 
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Assumption. For the expansions (p?<|), it will be assumed that, 

i>l j=0 

i>0 3=0 
Ni 

X±4> n = E E ^±^ J) (l - ^ W(l - r), 

i>0 j=0 

hold, where ' denotes differentiation with respect to p. 

The above calculus rules can be used together with the transport equations 
on ,y + to obtain equations for the coefficients <pn m the Ansatz (|38|). At 
each order (1 — r) p the solutions to these equations will be calculated. The 

(i j) 

conditions under which the coefficients (fy n can be pushed down to will then be 
investigated. The idea is then to compare the expansions obtained by means of 
this procedure with those arising from solving the transport equations on I. If 
an identification of the two different expansions is possible, then one can analyse 
the effect the regularity condition (|3^) has on the asymptotic expansions close 
to spatial infinity J? + . Now, if at each order (1 — r) p it is assumed that no 
log fe (l — t) terms are present in the lower terms, then it is possible to implement 
an inductive argument. The results of our forthcoming argument are summarised 
in the following theorem. 

Theorem 2. Let, 

i) the components of the spin 2 zero-rest-mass field 4> a bcd have an evolution 
expansion of the form 

Ni 

i>0 j=0 

n = 0, . . . , 4 with the coefficients <f>n at least C 1 in {r = 1, < p < po} for 
po a non-negative constant; 
ii) the calculus rules hold; 

Hi) (f>^ 0) (the radiation field) be C°° in {t = 1, < p < p a } c/+U 1+ ; 
iv) (f>i'°^ and (f)^ be bounded in {r = 1,0 < p < po}; 

v) the initial data </f> n | r= o, n = 0, ... ,4, be obtained by linearisation around 
flat space of asymptotically flat time symmetric initial data satisfying the 
Einstein vacuum constraints, and analytic in a neighbourhood of i; 

vi) furthermore, the initial data satisfy the regularity condition 

D(a s b s ■ ■ ■ D aibl b abcd) (i) = s = 0, 1, .. .; 



(39) 
(40) 
(41) 
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vii) be continuous in {r = 1,0 < p < po}; = 0, 

tial smoothness). 

Then the expansions of <f> n are in fact logarithm- free, that is, 

0„ = ^e O) (l-r) fe , 



where the coefficients cfin' ^ are C°° in {r = 1,0 < p < po}- 



, i + 2 f tangen- 



The base step (1 — r)°. In agreement with the Ansatz (38) let the leading 
terms of the spin 2 zero-rest-mass field be of the form, 

No Ni 

<t>n = £ ^ ln'(l - r) + £ £ <^')(1 - r) 1 - r), 

j=0 i>l j=0 

with n = . . . 4. Direct substitution into equations (||) show that, 

^ = ^ = ^ = 4 0J) = o, 

for j > 1. The coefficients (f>Q°\ 02°'°^ can be determined from the equa- 

tions (H) n = 0, 1, 2. The coefficients satisfy, 



, (0,0)/ 



,(0,0) _ , (0,0) _ 



0, 



(42) 
(43) 
(44) 



which can in principle be solved if the coefficient </> °'°' ) = </>o|,y+ (the radiation 
field) is given. It will be assumed that this is the case, and moreover, that it is 

C°°. Being ^q 0,0 " 1 smooth, one has that 



— J0,kr — 2—1 I 2—, 2—1 r 0,k;q,m 1 2q q-2 J P 

k>2 k>2 \q=2 m=0 / 



k 2q 



(45) 



with m complex constants. Using this expansion and equations ((42"l)-(ff4l) 

plus the further requirement of </>i and 03 to be smooth at I + one finds that, 

& = 0(A 0a = O(l), fo = 0{p). (46) 

The substitution of the Ansatz for the time development in equation (|^) n = 3, 
yields as a result the following hierarchy: 



6 (0,iVo)/ = ^ 



,(0,JV -1)/ 



-24°' 0) - ^ + p^' + X_4°' 0) = 0. 
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It can be solved from top to bottom, yielding: 



AO,N ) 

x(0,JVo- 
°4 



V 



6 (o,o) 



A ! 



X-^p~ 3 dp (47) 



where the hatted quantities denote functions which appear during the integration 
and thus, contain only angular dependence. The coefficients <p£ can be then 
"pushed down" to / + , but not in general in such a way that they are smooth at 
I + . The presence of terms of the form ln fc (l — r) in the Ansatz for the evolution 
yields as a result the presence of In °~ p terms in the corresponding expansions. 
This is an indication that if In p terms are present in the initial data then these 
will propagate in the evolution via some particular ln fc (l — r) terms. Logarithmic 
terms of the form ln J p do not appear in the normal expansions ( 34 ) . Thus, in 
order to be able of performing an identification of the expansions on and 
I some degree of tangential smoothness at I + will be assumed. A close look to 



formulae 



shows that one needs 



,(0,0) 



to have two continuous p-derivatives 



at I + . We will denote this by y A 

that $4° = for j > 2, and thus also 

a(0>! 



'" G Cp(I + ). This requirement now implies 



for j > 2. Note that the 
"integration constant" <t)f'^' does not need to vanish as the p 2 hip term in 4 ' ^ 
may be canceled with a similar one coming from p 2 j X^<f)f ^ p ^dp. Indeed, 

P~ 



the integral p 2 J X-cj)^' * 1 p 3 dp produces a term of the form p 2 In p from the p 2 



term in 



do,o) 



. As discussed previously, 



4°' 0) 



k>l 



(0) k 
3,kP ) 



(48) 



where the coefficients f^°l are calculated from the f^l of equation (|45|) by means 

of equations (fj^-([44"|). Setting (j)^' 1 ^ = X-f^\ i n accordance with our previous 
discussion one finds that, 



f (0) 



t 1} 



- X f (0) ^ 



5fV 



Y f(°) 

E ^~J3,k k 

k>l 
k^2 



(49) 
(50) 



One expects to be able to deduce the required degree of tangential smoothness 
form the smoothness properties of the initial data. 

Finally, one is now in the position of performing an identification of the ex- 
pansions obtained from the analysis of the transport equations on the cylinder, 
and those obtained from the transport equations at J? + . As discussed, both 
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expansions contain ln(l — r) terms. However, the ln(l — t)'s in the normal ex- 
pansions appear only at the highest spherical sector at each order. So do those 
in our asymptotic expansions: from equation ( |4g| ) one has, 

2 2q 

f(0) \^ \^ p(0) rp m 

J 3,2 — 2-^1 2-^1 r 3,2;q,m 1 2q q+\- 
q—1 m— 



And hence, 



2 2q 

A -/3,2 — ~2-^ t 'i,2;q,mP2q,q+2±2q q+2> 

q—1 m—0 



with (32q. q +2 = + 2)(cjr — 1). Thus, X_ cancels the sector given by q = 1, 
and as a consequence the only non-zero sector in (|9|) is g = 2. It is noticed 
by passing that X-4>^^ — whence it follows that (tangential smoothness 
assumed) X^cj)^ is bounded at I + . Hence, the leading term of the evolution 
expansion of 04 has the same form of the expansion given in ( |35| ) and (|36|). 
Consequently, if b a bcd(i) = (regularity condition at order s = 0) holds, we have 
shown that, 

Ni 

<t>n = +EE - r)* In'(l - r), (51) 

»>1 j=0 

n = 0, . ..,4. That is, there are no logarithms in the leading term of the expan- 
sions. 



The step (1 — r) 1 . The base induction step is somehow exceptional. In order 
to understand better the situation for the general case it is convenient to take a 
look to what happens at the order (1 — r) 1 . To begin with, let us assume that 
the analysis of base step (1 — r)° has been carried out, and no ln(l — r) are 
present at that order. Consequently, the expansion ( |5l| ) holds. 

Using equations ([?]) n — 0, 1, 2 and a similar approach to the one used in the 
base step, one readily finds that: 



0. 



for j > 1. This is a direct consequence of the fact that there are no logarithms 
at order (1 — r)°. Furthermore, 



6 (i,o) 



l r 

2 

1 - 

2 . 

1 " 

2 . 

l r 
2 



24°' 0) 



6 (o,o) 



4 (o,o)/ . 



-4°' 0) 



,a(o>°)' 

m 

6 (o,o), 4 



X 4 



j,(0,0) 

-01 

,(0,0) 



x 4 



4 o,o) 



(52) 
(53) 
(54) 
(55) 
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Thus, one can calculate the coefficients 



J. 1 ^ 



and 



directly 



from a knowledge of the coefficients at order p = 0. Using ( |46| ) in equations 
(p2|)-(^5|) one concludes that, 



(i.o) 



E &p k - 

k>\2-n\ 



The coefficients \ given in terms of the coefficients f!fl of the previous order. 
As in the base step, the logarithmic dependence comes from equation (||) with 
n = 3. From it one obtains, 



-30. 



(l,JVi-l) 



b<£> Nl) ' = o, 



,(l,iVl-l)/ 



= 0, 



-34 1,0) 



which again is solved from top to bottom so that, 



^-1)^1,^-1)^3 



r(i,iVi) 3, 



6 (i,o) 



1 

ivTI' 



X 1 ,^) 3 1 Ni 

>\ 'p°\n \ 



X- 



p dp. 



Whence one observes one more time the appearance of ln fe p dependence asso- 
ciated with the ln^ 1- ^! — r) in the Ansatz for the evolution. The expression 
p 3 J I_i^3 'p" 4 (i/) can only account for the cancellation of the p 3 lnp term in 



,(i,o) 



• j via the p 6 term in 4^'^- The coefficient 4>^' itself contains no logarith- 
mic dependence, as a consequence of equation ( |55j ) and the smoothness of the 
coefficients calculated at the base step p = 0. Following the spirit of previous dis- 
cussions, we require (f>^' ' 6 C'p{I + ). Therefore $4 = for j = 2, . . . , N\, and 

furthermore = for j = 2, . . . , N x . So, one is left with only two non-zero 

coefficients. Namely, 



(1,0) 



A_ 



dp 



Now, writing 



*!>i 



(56) 
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one has to set 



h - x -h,3' 



(57) 



in order to eliminate the remaining logarithmic term in V4 

In a similar way to what happened in the base step, one can further show 
that the logarithmic dependence is only found at the highest spherical harmonic 
sector. However, in this case and also in the general case, the analysis is a bit 
more elaborate. From equations (E9I) and (pH) one finds that: 



6 (i,0) _ 1 
"3 — 2 



(0) 

2 P 



k>l 



k>3 



The spin weight of ^3 is —1, and accordingly, 

fc 2q 

AO) _ p(°) r m 

J3,k ~ Z^ 3,k-q, m - L 2q q+1- 
q=l m=0 

Thus, at the end of the day one has the following rather complicate- looking 
expression: 

4 1,0) = £/iV 

fc>2 



fc 2q 



£ E E ( is± » 11 - ( fc + 1) 



fc>l q=l m=0 
fc^2 

The term (q + 2)(q — 1) — (fc + l)(fc — 2) vanishes whenever q = fc — 1. Hence, if 
one writes 

fc 2q 
\^ \^ p (l) T m 
■>3,fc ~~ Z^ Z^ 3,k;q,m^2q q+1) 
q—1 m—0 

for fc = 1, 2, . . . . Therefore, the term in X. 



one has F, 



(1) 



3,fc:fc — l.m 



contains only q — 3 and g = 1 spherical harmonic dependence, i.e. 

2 6 

•>3,3 — Z^ 3,3;l,m- t 2 2 ' Z^ 3,3;3,m J 6 
771=0 m=0 

Therefore, using (p7p one finally arrives to: 



l f> 4- 



4 M) 



= - V E ^UmT, 



6 5- 



(58) 



m=0 
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Again, the expansions have the same form of ( p5| ) and (|3^) with s* = 0. If the reg- 
ularity condition is satisfied up to s = 1, i.e. if fr a &c<i(*) = and D^ aibl b abcd ^ (i) = 
then 

Ni 

K = + <^>°>(1 - r) + " r)' W(l - r). 

The general step (1 — t) p + 1 . The procedures of the general step are lengthier 
but nevertheless of a similar nature to those of the analysis of the (1 — r) 1 terms. 
We begin by assuming that a similar analysis to that one carried for (1 — r) 1 has 
been carried up to (1 — r) p inclusive. Accordingly, the components of the field 
are assumed to have an expansion of the form, 

p N k 
0„ = $>^ O) (l-r) fc + £ £^>(1-t)*1ii'(1-t). (59) 

k=0 k>p+l j=0 

Ik 0) 

The coefficients (jr n ' with k < p are C°°, for they have been constructed out 
of the radiation field </>o°'°' ) in a way that preserves the smoothness. Hence, 

fe>|2— n\ 

where the coefficient f n ,k contains all the angular dependence. In particular, for 
the component </>3 they are such that, 

k 2q 

f(p) _ \ " \ jp(j>) rp m 

J3,k — 2-^i 2-^i 3,k;q,m ± 2q q+V 
q—1 m—0 

with coefficients f!: p ) „ such that 

0,^,0, Tit 

jpip) _ jp(p) _ _ P (p) _ n ( RCl \ 

r 3,k;b-l,m ~ r 3,k;k-2,m — ' ' ' — r 3 1 k-k~p,m ~ u ' 

as suggested by the analysis of the (1 — r)° and (1 — r) 1 cases. Substitution of 
the Ansatz (|5^) in equations (Q) with n = 0,1,2,3 yields a hierarchy of N p+ i + 1 
equations. A similar analysis to that performed at order (1 — r) readily shows 
that 

^+^"> = <^ +i ^ = = = o, 

for 1 <j < Np+i. One also obtains a set of 4 recurrence relations which allows 
to calculate the coefficients </>n P+1,0 ' ) n = 0, 1, 2, 3 from (in principle known) lower 
order terms. In future discussions only the expression for </>g P+1 ' ^ will be needed, 



,(p+l,0) _ 1 
'3 - 2(p+l) 



( P -l)# 0) -^- 0) '+A + 4^1. (61) 
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The coefficients (jif^ 1 ^ are on the other hand calculated using equation (|J) with 
n = 3. Similar computations to those described in the (1 — r)° and (1 — r) 1 steps 
lead to: 

,(p+l,N p+1 ) _ 2iP+hN P+ i) „ +3 
— Pi 



^+1,0) = _t_ftp+l,N p+l)pP+3 lnNp+1 p+ ... + ^+l,0) pP+ 3 

- ff+ 3 y X^^p-^dp. (62) 

Again, the term p p+3 J X_(j)^ +1 ' ^ p~^~ v dp in the hierarchy ( |62| ) can be only 
used to cancel out the cpf pP+ 3 hip term in the expression for the coeffi- 
cient 4 P+1 ' 0) - Thus, in order to perform an identification at I + , some tangen- 
tial smoothness will be required. More precisely, we will demand </>4 P+1,0 ' ) G 
CP +3 (I+). This assumption implies that $f +1 ' Np+l) = ... = 4>{ p+1 ' 2) = 0, with 

the further consequence of yielding 0^ P+1 ' A ' P+1 ' ) = . . . = <^| +1 ' 2 ' = 0. Hence as 
before, we are left with only two non-zero coefficients: 

^+M) = ^ + i,iy+ 3) (63) 
4p +1 >°) = ^ +1 <V +3 Inp + ^ +1 'V +3 

- p p+3 [ X_$ +1 ' Q) p~ 4 ~rdp. (64) 



The coefficient 0g P+1 '°' ) is smooth, as it is constructed from the </>L P '°^ 's using 



equation (|61j) . Consequently, 

E-& + V, (65) 



5 (p- no; = /•'/' - J i 
fc>i 



with the coefficients f^ 1 ^ containing only angular dependence. Hence, one sees 
that: 

pr>+ 3 ( \_^t m p-^dp = X_f^l f 7+ 3 \np+ J2 



fe>i 

fc#p+3 



In accordance with our previous discussion one sets: 

;Kp+i,i) _ v f(P+l) 
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In this way we have accounted for the remaining logarithm in cj) l f +1 ' ^ — see 
equation (|6^). Therefore, at the end of the day one has that: 

,(p+l,l) _ Y Ap+l) p+3 
'Pi — jX -J3,p+3P ' 



(p+1,0) _ -7O+1.0) p+3 x -fj P k 1) n k 

k>2 



k^p+3 

The component fa is the most irregular one of the field, and thus, the ln(l — 
t) terms will appear firstly there. Carrying out the previous discussion to the 
following orders, one finds that the first ln(l — r) terms in fa, fa, fa and fa 
appear at orders (1 — r) p+2 , (1 — t) p+3 , (1 — r) p+4 and (1 — r) p+5 respectively. 

Finally, in order to make use of the regularity condition ([33]) it is again 
necessary to identify the ln(l — r)'s appearing in our expansions with those 
coming from the analysis of the transport equations on the cylinder I. Using 
equations (^) and (|5|) one finds 



<, (p+i,o) 



P+l 



fw^+Efo-i-*)^- E n£#p* 

k>l k>l 
k^p+2 



,(p+l,0) _ 1 
5 3 - p+l 



One can make explicit the spherical harmonics dependence: 

+ E EE i^^ 1 + (p - * - 1)) 

fc>l <j=l m=0 

where the coefficients Fg m satisfy (|60|). The expression (fc — p + l)(fe — p — 
2) — (q + 2)(g — 1) is zero whenever q = p — k and/or q — k—p — 1. Hence, writing 

k 2q 

f (P+l) \^ p(P +1 ) T m 

J3,fe ~~ 3.k;q,m- L 2q g+D 

q— 1 m— 



then one finds that 



p(p+i) _ n 

-^fcjfc-p-l.m — u ' 



Furthermore, using the induction hypothesis ( j60| ) in equation ( p6[ ) one finds that 

p(P+l) _ p(P+l) = p(P+l) p(P+l) _ Q 

3,k;k — l,m 3,fc;fc — 2,m ' ' ' 3,fc;fc— p,m k;k— p— l,m 

In particular if fc = p + 3 then X_/p + 1^p p+3 m P wu l only contain spherical 
harmonic dependence at the q = p + 3 sector. The other possibly remaining 
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sector, the q — 1 one, is annihilated by the AT_ derivative. Accordingly, it has 
been proved that 

2p+6 

,(p+l,l) _ _ R p+3 P (p+1) ^ rn 

< t ) 4 — l J 2p+6,p+5P 2_^i r 3,p+3;p+3, m 1 2p+6 p+5' 

Hence, the expansions we have so far obtained are similar to those given in ([55]) 
and (|3^) . Finally, if the regularity condition (|3^) holds up to order p+1 then 
one has, 

P+1 N k 

*n = £^ 0, (i-r)*+ E E^a-^a-r). 

fc=0 fc>p+2 j=0 

Concluding remarks. Some final remarks come now into place. One would 
expect to be able to deduce the smoothness of 0o on null infinity, the bound- 
edness of <f>\ and 03 and the tangential smoothness of </>4 from the smoothness 
of the initial data. This would require, in principle, the implementation of some 
energy estimates which would allow us to reach null infinity (see the discus- 
sion in §3). This is at the time of writing still an open problem. The present 
analysis is nevertheless valuable in the sense that it focuses our attention on 
the facts/properties one should be able to prove, and their interconnections. It 
should be, in principle, possible to undertake an analogue of the above discussion 
for the non-linear gravitational field. The analysis of the linearised gravitational 
field has shown us the way to lead. These matters will be the subject of future 
work. 
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